use ipums, clear

// Muslims as share of total
graph bar muslim algeria morocco tunisia turkey, bargap(20) bar(1, fcolor(blue) lcolor(blue)) ///
	bar(2, fcolor(blue) lcolor(blue)) bar(3, fcolor(blue) lcolor(blue)) bar(4, fcolor(blue) lcolor(blue)) bar(5, fcolor(blue) lcolor(blue)) ///
	intensity(*0.7) lintensity(*1) showyvars legend(off) title("Father's birthplace") ///
	ysca(titlegap(2)) ylabel(, nogrid labsize(small)) ytitle("Share of native-born women aged 20-33", size(small)) ///
	yvaroptions(relabel(1 "Muslim-majority" 2 "Algeria" 3 "Morocco" 4 "Tunisia" 5 "Turkey") label(labsize(small)))

	
// Comparison of education and LFP, Muslim vs non-Muslim
collapse secondary university lfp employed, by(muslim)

tsset muslim
foreach x in secondary university lfp employed {
	gen D`x'=D.`x'
}

graph bar Dsecondary Duniversity Dlfp Demployed, bargap(20) bar(1, fcolor(red) lcolor(red)) ///
	bar(2, fcolor(red) lcolor(red)) bar(3, fcolor(red) lcolor(red)) bar(4, fcolor(red) lcolor(red)) bar(5, fcolor(red) lcolor(red)) ///
	intensity(*0.7) lintensity(*1) showyvars legend(off) title("Education and labor force participation") yline(0, lcolor(gs10)) ///
	ysca(titlegap(2)) ylabel(, nogrid labsize(small)) ytitle("Difference Muslim - Non-Muslim origin women aged 20-33", size(small)) ///
	yvaroptions(relabel(1 "Completed secondary" 2 "University degree" 3 "In labor force" 4 "Employed") label(labsize(small)))
